00 
ON 
ON 



ON 



X 



Strange hyperon and antihyperon production from 
quark and string-rope matter 

P. Csizmadia 1 , P. Levai 1,2 t S.E. Vance 2 , 
T.S. Biro 1 , M. Gyulassy 2 , J. Zimanyi 1 



1 KFKI Research Institute for Particle and Nuclear Physics, 
Qh 1 P- O. Box 49, Budapest, 1525, Hungary 

2 Department of Physics, Columbia University, 
New York, NY, 10027, USA 

18 September 1998 



> 
in 

ON . 

Abstract 

00 

Hyperon and antihyperon production is investigated using two microscopical 
,1^ models: (1) the fast hadronization of quark matter as given by the ALCOR model; 

(2) string formation and fragmentation as in the HIJING/B model. We calculate 
the particle numbers and momentum distributions for Pb+Pb collisions at CERN 
SPS energies in order to compare the two models with each other and with the 
available experimental data. We show that these two theoretical approaches give 
similar yields for the hyperons, but strongly differ for antihyperons. 



1 Introduction 

The quest to produce the Quark-Gluon Plasma (QGP) at the CERN SPS has reached its 
final stage with the Pb+Pb run at Eb eam = 158 AGeV. Since strangeness enhancement 
in heavy ion collisions relative to pp collision has been predicted as a QGP signal |T], [2j , 
the measurement of strange and antistrange particles, especially the hyperons and the 
antihyperons, has received much attention (see Refs. PH-RIll). Several thermal models 
were invented to reproduce earlier experimental data at the AGS and the SPS energies |12 



13| , IHH . However, only extended versions which introduce more parameters are capable of 



fitting the new experimental data which particularly include the hyperon and antihyperon 
yields |15| . Although these thermal models are capable of describing the gross yields and 
thermodynamic properties which are observed in the final state of the collisions, they are 
unable to provide a description of the various underlying microscopical mechanisms which 
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may be involved in the initial stages following the collision of the two nuclei. Thus, to test 
for the formation of a QGP or any other possible initial phase, we choose two microscopical 
models which describe the evolution of heavy ion collisions from two different viewpoints; 
a plasma vs. a hadronic scenario. 

The first model assumes the formation of a deconfined state (namely the Constituent 
Quark Plasma (CQP) [0, 17]) and then describes its hadronization. The applied ALge- 



braic Coalescence Rehadronization (ALCOR) model [ff8L O and Microscopical Coales- 



cence Rehadronization (MICOR) model p0[ (as well the time-dependent Transchemistry 
Model jn| [T7|]) are all based on this assumption. 



The second model begins with a description of hadron-hadron collisions using string 
phenomenology. The strings are then fragmented to produce the final observed hadron 



spectrum. The HIJING /B [^TJ event generator (an extension of HIJING |22] event gen- 



erator to include a baryon number transport mechanism) is one version of this model. 
The Relativistic Quantum Molecular Dynamics (RQMD) model is another variant of this 
model p3|. However, we note that HIJING and HIJING/B greatly differ from RQMD in 



that they include minijets. 

Although both of these models are able to reproduce pion, kaon and proton data, they 
differ in their predictions of the hyperons and the antihyperons. Thus the underlying 
assumptions of these models can be investigated by comparing their theoretical predictions 
for the hyperon and the antihyperon production with the experimental data. 

Before discussing our calculations, we argue that secondary hadronic collisions neg- 
ligibly effect the measured hadronic spectra at SPS energies. We note that secondary 
hadronic interactions are needed to describe the lower energy heavy ion collisions (at SIS 
and AGS) as seen in the successful use of hadronic cascade models in reproducing data at 
these energies. However, at Eb eam « 100-200 GeV/nucleon we expect that the collisions 
involve the quarks of the nucleon. This introduces uncertain hadronic formation time with 
the production of a hadron from the intermediate quark matter or string/rope resonance. 
This time is on the order of the life-time of the early highly excited system. This idea is 
just the old argument given by Pomeranchuk and Feinberg p4| who claimed that particles 
only really materialize, not at the early collision stage, but at low energy densities where 
they freeze-out and become interactionless; the "free separation" of particles. 

The HIJING and HIJING/B models use this assumption, since fragmentation occurs 
after the nuclei have finished colliding. (We will also display the 'no-scattering' results of 
RQMD from Ref. 0.) 

The ALCOR and MICOR models assume that the hadronization of deconfined matter 
yields an excited hadronic gas which includes excited resonances, namely the vector meson 
nonet and the baryon decuplet beyond the lowest lying pseudoscalar nonet and baryon 
octet. The resonances decay into the stable hadrons. All of these models, except RQMD 
neglect final state interactions. 

Since the final hadron production of these models is directly related to the underly- 
ing microscopical mechanisms, we compare their predictions of hyperon and antihyperon 
formation with each other and with the available experimental data. 
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2 The ALCOR and MICOR models 

In Pb+Pb collisions at the CERN SPS energies, it is assumed that the incoming nucleons 
disintegrate into a form of deconfmed matter. Originally the appearance of the perturba- 
tive QGP phase was expected. However, around the critical temperature, T c « 150 — 200 
MeV, perturbative QCD is unable to describe the characteristic processes completely, 
thus non-perturbative or phenomenological descriptions should be introduced for the de- 
confined matter. The non-perturbative descriptions can be motivated by lattice QCD 
calculations of the equation of state near critical temperature T c . Recent analysis indi- 
cates p5j that deconfmed matter in its equilibrium state can be described by massive 



quarks and gluons. These gluon-like quasi-particles are heavier than the quarks (thus 
their number is suppressed if the system is thermalized) and the effective quark masses 
are close to the constituent quark masses of the additive quark model We call this 



deconfmed state the Constituent Quark Plasma (CQP). The CQP state is dominated by 
quarks and antiquarks, where some combine to form colorless hadron-like objects. The 
presence of such prehadrons in the deconfmed phase is confirmed by lattice-QCD calcu- 
lations of hadronic correlation functions [^7j. It is assumed that these prehadrons can 
escape from the strongly interacting region to become mostly excited "real" hadrons. 

The ALCOR model describes this microscopical hadronization process using quark 
coalescence (for a detailed description see Refs.J18|, The coalescence rates depend on 
the temperature T, on the thermal quark masses M q and M s , and on the strength of the 
coalescence, which can be characterized by an effective a e ff coupling in a Coulomb-like 
quark-quark interaction potential. From the point of view of hyperon and antihyperon 
production it is important to note, that the applied coalescence cross section is pro- 
portional to the mass of produced prehadron, <J l C oi[^ h oc m\ (see Ref. []18|). Thus the 
production of heavier particles, especially the hyperons (antihyperons) containing the 
heavy strange quark (antiquark) is relatively enhanced. To describe the hadron forma- 
tion completely, two more parameters are needed: the number of newly produced light 
quarks N qq = N U u = and number of strange quarks N S g. The values of N qq and N s % 
are scaled by A 13 , where (3 ~ 1. These input parameters, combined with some initial 
geometry and a longitudinal Bjorken scaling are the basis of the ALCOR model which 
calculates the total number of produced hadrons for different heavy ion collisions. 

In the ALCOR model the coalescence rates are determined from multidimensional 
integrals which can be calculated exactly for the case of Bjorken scaling. To obtain the 
differential momentum spectra of the produced hadrons, integrals need to be separated 
and differential coalescence rates need to be evaluated numerically. The introduction of 
numerical integrals allows us to consider a more complicated flow pattern, in the present 
case a constant transverse flow characterized by v tr , and a finite longitudinal extension 
characterized by the maximal space-time rapidity f] ma x- A new version of the ALCOR 



model, the Microscopical Coalescence Rehadronization (MICOR) model pOfl , includes this 
new flow pattern and calculate the differential momentum spectra. Here, the microscopi- 
cal quark coalescence produces off-shell prehadrons and we assume that these prehadrons 
escape the deconfmed region conserving their velocity to become real on-shell hadrons. 
(Because of the small deviation between prehadron and real hadron masses, this assump- 
tion works well.) After calculating the differential rates, the absolute norms of the particle 
production are determined by a set of equations based on quark and antiquark number 



3 



conservation, similarly to ALCOR. The MICOR model is able to determine the differen- 
tial hadron spectra which arises from the hadronization of the CQP state. As in ALCOR, 
final state hadronic rescattering is neglected. 

The advantage of using the ALCOR and MICOR models is that the phase transition 
can be followed microscopically. In this way a phenomenological quark-quark interaction 
determines the formation of a chemically out-of equilibrium hadronic gas. 



3 The HIJING/B model 



As has been summarized elsewhere [|2| , the HIJING Monte Carlo event generator models 
hadronic interactions by combining low px multistring phenomenology with the pertur- 
bative QCD (pQCD) processes. The HIJING/B event generator extends the HIJING 
model by introducing a novel, non-perturbative gluon junction mechanism in order to 
study baryon number transport in heavy-ion collisions. As shown below, this novel non- 
perturbative junction mechanism strongly influences hyperon production. 

This gluonic mechanism is motivated from the non-perturbative gluon field config- 
uration (the baryon junction) that appears when writing the simplest gauge invariant 
operator for the baryon in SU C (3); 



B = e ijk 



Pexp(igJ dx^A^j q{x x ) 
Pexp(igJ dx^A^j q(x 2 ) P exp (ig J dx^A^j q(x 3 



(1) 



Here, the baryon junction is the vertex at xj where the three gluon Wilson lines link the 
three valence quarks to form the gauge invariant non-local operator. In a highly excited 
baryonic state, the Wilson lines represent color flux tubes. When these strings fragment 
via qq production, the resulting baryon will be composed of the three sea quarks which 
are linked to the junction while the original valence quarks will emerge as constituents 
of three leading mesons. Being a gluonic configuration, it was proposed [p8] that the 



junction could be easily transported into the mid-rapidity region in hadronic interactions. 

In Regge phenomenology, allowing for the possibility of baryon junctions in scattering 
is taken into account by introducing additional Regge trajectories (M J ). As was shown 
by Kharzeev [28fl, the exchange of these trajectories in p + p interactions leads to events 



where one of the baryon junctions is stopped in the mid-rapidity region. Once the baryon 
junction is stopped, the three valence quarks continue forming three beam jets. These 
events have a different energy, rapidity and multiplicity dependence compared to the 
normal Pomeron two string event. In particular, the produced baryon has a cosh(?//2) 
rapidity dependence and a l/tfs energy dependence. In addition, the new string config- 
uration leads to a multiplicity enhancement of a factor of 5/4 (counting the two + three 
jets as compared with the two + two jets). Since the baryon which is resolved around the 
junction in the mid-rapidity region is composed of three sea quarks, hyperon production 
and the < pi > of the final baryon are enhanced by factors of 3. The enhancement factor 
of 3 of the strangeness content of the baryon allows for the unique possibility of producing 
Q~ (S = —3) baryons. 
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In HIJING/B PI], the baryon junction mechanism is implemented using a "Y" string 
configuration for the excited baryon with an effective cross section of ~ 9 mb. The above 
junction dynamics does not provide a mechanism for antibaryons production and further 
studies which address baryon pair production are underway. In their present version, 
neither HIJING nor HIJING/B include final state interactions. Thus, in one variant of 



HIJING/B we simulate multiple final state interactions using the concept of "ropes" [29 
increasing the energy per unit length of the strings from 1.0 GeV/fm to 1.45 GeV/fm in the 
heavy nuclear collisions. Increasing the effective string tension increases the strangeness 
production and the < px >■ 



4 Results on PbPb collision 

In the ALCOR model, we recalculate the hadron production as compared with earlier 
publications |1| [19]]. In this calculation, we considered the feeding of the detected As 
from S and Q decay channels and the undistinguishable S°. We also introduce the A — 
like = A + S° + H~ + H° + Q~ species, which was not done earlier. This number will be 
compared to the experimentally measured number of "A" -particles, as was done similarly 
with "A°-like" particles. In this case the input parameters for newly produced quark pairs 
are N qq = 391 and N S g = 172 (which means g s = N S ^/(N U ^ + N^) = 0.22) and for the 
effective quark-quark interaction, a e // = 0.97. The results are displayed in Table 1. 

The particle yields were also calculated in HIJING/B with ropes where the energy 
per unit length of the string was increased to k = 1.45 GeV/fm. The total numbers are 



shown in Table 1 along with the RQMD results [[23 1. 

When comparing the two calculations, ALCOR and HIJING/B predict approxi- 
mately the same numbers of mesons, baryons and non-strange antibaryons; approximately 
reproducing the experimental data. However, significant differences can be seen in the 
predictions of strange antibaryons and the Q. These differences are directly related to the 
assumed microscopical hadronization mechanisms and thus reveal the applicability of the 
different assumptions. These differences can be investigated via particle ratios, which are 
displayed in Table 2 and 3. In Table 2 and in Table 3, the experimental results obtained 
from NA49 Collaboration and the WA97 Collaboration are shown, respectively. We note 
that these measured ratios are valid in certain rapidity windows, and thus only provide 
an approximate comparison. 

We also investigate the energy dependence of our results. Particle abundances were 
calculated at 100, 200 and 300 AGeV for Pb+Pb collisions using HIJING/B with ropes 
and the ALCOR model. In the ALCOR model we kept the g s and a e ff parameters at 
their value obtained by fitting the data at 158 AGeV, and only changed the number of 
newly produced light quarks, N qq = 354, 442, 496, respectively, to reproduce the h~ num- 
ber as obtained in HIJING/B. The results of the energy dependence of strange particle 
ratios from these calculations are shown in Figure 1. In addition, we show the experi- 
mental results from the NA49 and WA97 experiments. Although both models reveal a 
slight energy dependence, only the ALCOR model was able to consistently reproduce the 
observed values. 
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Table 1: Total hadron multiplicities for Pb + Pb collision at bombarding energy 158 
GeV/nucleon. The displayed experimental results are obtained or estimated ({}) from 
the publications of NA49 Collaboration: a is from UlQ |; b is from [0; c is from [11]]; d is 



from ||; * is estimated from {K~} and {K$}. Theoretical results are calculated in the 
ALCOR and HIJING/B models and it was included the results of RQMD model ("ropes 
+ no re-scattering" version) [p3|| . We introduced A — like = A + S° + 2~ + 2° + Q~. 
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Table 2. Strange baryon and antibaryon ratios measured by NA49 Coll. for Pb ± Pb 
collision at 158 GeV/nucl bombarding energy in the momentum region 3.1 > y > 4.1, 
Pt > 0. We calculated the ratios of total numbers from Table 1. Here A (A) means 
"A°-like" ("A°-like") species. 
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Table 3. The ratios of strange baryons and antibaryons, measured by the WA97 Col- 
laboration |J for Pb ± Pb collision at 158 GeV/nucl bombarding energy at px > GeV. 
We calculated the ratios of total numbers from Table 1. Here A(A) means "A°-like" 
("AMike") species. 




Figure 1 . The measured ratios of baryons and antibaryons in the PbPb 
collision at 158 AGeV energy from the NA49 Coll. (stars) g, [Kj and 
the WA97 Coll. (squares) 0. The theoretical results from ALCOR 
(dots) and HIJING/B (triangles) are displayed at 100, 200 and 300 AGeV 
energy. 
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The transverse momentum distributions of the produced hadrons in the mid-rapidity 
region were also calculated using the MICOR model. For this calculation, two initial 
conditions were considered. Since the ALCOR model is not very sensitive to the initial 
temperature Tcqp, in the range 150 MeV < Tqqp < 180 MeV, two set of parameters were 
used in the MICOR calculation: I. T CQP = 150 MeV, v tr = 0.67; II. T CQP = 175 MeV, 
v tr = 0.57. In case I, the pion and proton spectra is reproduced at midrapidity. In case 
II, good agreement was found for the heavier hyperons. 



Baryons — Pb+Pb, y=0, Set II — Anti— baryons 




m T — m m T — m 



Figure 2. The transverse momentum distribution of proton, A, 2~ 
and their antiparticles in Pb+Pb collision in midrapidity, produced by 
the MICOR model with initial condition set II. The full curves are the 
results of calculation and the dashed lines indicate experimental slopes 
measured by NA49 Coll. @. 

In the MICOR model, the non-thermal momentum distributions of the emerging ex- 
cited hadrons can be calculated. These distributions can be characterized by a non- 
Boltzmann distribution, / oc (m/E) s ■ exp(—E/T), where 5 ~ 2. What is remarkable is 
that after resonance decay the obtained hadronic spectra can be very well fitted by pure 
exponential functions |20| , with the effective slopes closely matching the experimental val- 
ues. In Figure 2, we display the resulting momentum spectra using the values from case 
II for the proton, A, S~ and their antiparticles. Here, the solid lines show the MICOR 
results, which can be parametrized by the exponential fit: 

d JL— = C ■ e - mT ^ff . (2) 

mTdrriTdy 

These results can be compared to the dashed lines which represent the slopes as 
measured by the NA49 Collaboration 0, ||. Here, we show that after two completely 
non-equilibrium steps (quark coalescence and decay of excited hadrons) the final hadron 
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spectra is approximately exponential. Our result in Figure 2 demonstrates that non- 
equilibrium processes are able to produce particle spectra which can also be interpreted 
as thermal. We note that baryons and antibaryons have approximately the same slope, 
as has been observed. 

In Figure 3, we display the effective slopes T e // as calculated using the MICOR model 
(using cases I and II) and as measured by the NA49 @, g, NA44 @ and WA97 Collabora- 
tions H H . The possible interpretation of our MICOR results is the following: parameter 
sets I and II are so close to each other that a mutual hadronization temperature and 
transverse flow inbetween (namely T C qp ~ 170MeV and v tr ~ 0.6) could be assumed 
for the produced deconfined state as the dotted line shows in Figure 3 (see Ref. [[|). 
The deviation of the recent experimental data from this line indicates that hyperons (an- 
tihyperons) which contain mostly strange (antistrange) quarks escaped the light-quark 
dominated CQP state at a slightly earlier time, when the temperature is higher and the 
transverse flow is smaller. Particles containing mostly light quarks are formed in the 
later stage of the expanding CQP matter and thus their spectra indicate a lower tem- 
perature but higher transverse flow upon hadronization. We note that there is a similar, 
interacting hadronic explanation which was obtained from an RQMD-based analysis 
which assumes secondary interactions with different strengths for strange and non-strange 
particles. We hope that final data on Q~ and f2 + production will help us to make our 
conclusion more definite. 
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Figure 3. The effective slope parameters, T e ff, in the midrapidity region 
for Pb+Pb collision at 158 AGeV, obtained from the MICOR model and 
from the different experiments. The dotted line indicates the appearance 
of a common transverse flow in the hadronic phase, see Ref. @. 
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5 Conclusion 



We investigate hyperon and antihyperon production in Pb+Pb collision at CERN SPS 
energies using two different scenarios. The formation of a Constituent Quark Plasma and 
the modeling of its hadronization as in the ALCOR and MICOR models was used to 
calculate the particle numbers and transverse momentum distributions. String formation 
and fragmentation as in the HIJING/B with ropes model was also used to calculate 
particle yields. In the case of meson, baryon and non-strange antibaryon yields, both 
models predict similar particle numbers. However, only the CQP scenario could reproduce 
the antihyperon data. This can be attributed to the coalescence of heavier strange quarks 
which is favoured in the ALCOR model. The inability of HIJING/B to reproduce the 
antihyperon data reveals the need for a new hadronic mechanisms (similar to the baryon 
junction) for antihyperon production. 

In addition to reproducing the antihyperon yields, the results of the CQP-based MI- 
COR model point to a scenario where a mutual transverse flow was developed in the 
deconfined phase and that hadronization occured when Tcqp ~ 170 MeV and Vt r ~ 0.6. 
Here, hyperons and antihyperons may escape a little bit earlier from the deconfined re- 
gion, retaining a slightly larger temperature and a smaller collective transverse flow. The 
other hadrons are then formed shortly thereafter. 
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